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1. Introduction 



Direct detection of extrasolar planets by imaging is an area of major interest in the 
study of other solar systems, opening up wide-orbit non-transiting planets to photometric 
and spectroscopic characterization. Currently, only mass ive giant planes in wide orbits ca n 



be detected against the bri ght light of their parent star (IMarois et al 



2008 



2008a 



Kalas et al. 



Lagrange et al. 



20091 ) at contrast levels of approximately 10 6 and separations of 



one arcsecond or more. The ultimate goal is the direct detection of terre strial planets at 
contr ast levels of 10 10 and separations of a tenth of an arcsecond or less (jBeichman et al. 



19991 1 . The next major step will be dedicated grou nd-based high-contrast adaptive optics 



systems such as th e Gemini Planet Imager (GPI) ([Macintosh et al 



( IBeuzit et al 



2008|). 



20081 ) and SPHERE 



A crucial component of such a system is the coronagraph - broadly defined as a 
device intended to block the diffraction of the coherent part of the optical wavefront. 
Considerable effort has gone into the design of coronagraphs suitable for the extremely high 
contrast space application; somew hat less effort h a s gon e into optimizing ground-based AO 



coronagraphs, though recent work 



Martinez et al 



(120081 ) has appeared. In many ways, the 



latter is a more challenging problem. Ground-based telescopes almost inevitably have large 
secondary mirror obscurations. Another challenge is posed by the nature of the wavefront 
presented to a ground-based coronagraph. Even after adaptive optics correction, the 
wavefront will include a large amount of high-spatial-frequency error - tens of nanometers 
rms, both dynamic and static - just outside the "dark hole" frequency cutoff of the AO 
system. Due to the constant motion of the atmosphere, there will also be smaller (but 
still significant) amounts of dynamic wavefront error within the controlled range. Any 
ground-based coronagraph must be robust against these wavefront errors. These wavefront 
errors produce an instantaneous contrast floor of approximately 10 5 from 0.2-0.8 arcseconds, 
but the performance of the coronagraph must be much better than this so that "pinned 



2 



speckle" cross-talk between the small static com ponent of the wavefront error and the 



coronagraph residuals do not limit final contrast (lAime and Soummer 



2004 



Soummer et al. 



20071 ). Multiwavelength imaging speckle-rejection techniques can somewhat relax these 
requirements, but even in that case, cross-terms between the residual diffraction and the 
quasi- static chromatic and achromatic wavefront error can significantly limit final contrast, 
particularly since the residual diffraction will often have a very different chromaticity than 
the speckles induced by atmospheric wavefront error. 

The baseline coronagraph design for GPI is an apodized pupil Lyot coronagraph 
(APLC). An APLC is structurally similar to the Lyot coronagraph: following Fig. [TJ it 
starts with an aperture at a pupil plane (plane 1), a hard-edged occulting mask at an image 
plane (plane 2) and a second pupil plane with a Lyot stop (plane 3), and finally takes 
images at the second image plane (plane 4). Unlike the Lyot coronagraph, however, the 
aperture is not open but apodized; the apodization is a generalized prolate spheroidal wave 
function, chosen together with the focal plane mask and t he L yot stop to optim ally remove 



Soummer et al. 



(|2003j) and 



Soummerl (120051 ) for further 



on-axis light from the system. (See 
information.) 

The apodized pupil is shown in Fig. [2j Several methods to produce this apodization 
have been investiga ted and preliminary results have been obtained with 2/xm chromium 



microdots on glass (IThomas et al. 



20081 ). These dots are placed in a random pattern to 



attenuate the light passing through the pupil by the amount specified by th e apodization. 



The current performance of the apodizer has not been as good as expected (IThomas et al. 



20081 ); it is suspected some of this is due to the small size of the dots, which are on the same 
order as the wavelength. For example, light passing through a grid of conducting holes with 
size similar to the wavelength will experience a phase delay, which results in a apodization- 
and wavelength-dependent phase error across the pupil. A differ ent randomization m ethod 



with larger dots has been investigated to mitigate these effects ( IMartinez et al 



20091 ). An 
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alternative approach to mitigate the effects of small microdots is to use shaped pupils in 
the first pupil plane of the APLC; we present simulations of this approach in Section [3J 



2. Theory 



The electric field at the image plane following a normally-incident plane wave on an 
apodized, radially-symmetric pupil with transmission function A(r) may be written: 



E im {p) 



2tt 



R 



E A{r)J (^f ) 
Jo \ A J 



T 



Shaped pupils are masks placed in the pupil plane of a coronagraph, whose apertures 
have been designed, usually through optimization techniques, to produce r egions of high 



contrast in the following image plane. They com e in a variety of designs (IKasdin et al. 



2003b 



Vanderbei et al. 



2003 



Kasdin et al. 



2003al ) , but for the purposes of GPI, we wish to 



consi der one type in particular, the starshaped mask. A starshaped mask (j Vanderbei et al. 



20031 ) can be thoug 



it of as a binary appro ximation to an azimuthally-symmetric apodized 



pupil. As shown in 



Vanderbei et al. 



(120031 ). the electric field in the image plane following a 
starshaped mask can be written as the electric field from the apodized pupil plus a series of 
perturbation terms: 



im, binary 



(P,<f>) = E im (p) 



+ 
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x (2cos(jW(0- tt/2))) 



(2) 



where A(r) is the apodization profile on [0, 1] with 1 letting all light through, and blocking 
all light, and N is the number of holes, r here is the focal plane coordinate, (p, 0) are the 
image plane coordinates, and R is the radius of the pupil. The matching is accomplished 
by noting that a circle drawn with radius r, centered at the center of the mask, would have 
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a certain fraction of its circumference open, as it passes over holes in the mask. If this 
fraction is set to be A(r), the above equation results. (See Fig. [3] for an image of the holes.) 

The advantage of starshaped masks comes from the fact that the perturbation terms 
can be made very small for small p when the number of holes is large. This has two 
effects — first, it means that the outer working angle (OWA) of the shaped pupil can be 
pushed outward by expanding the number of holes, with little effect on the field at smaller 
angles, and second, it means that an image plane mask (such as that in the APLC) blocks 
very little of the electric field from the perturbation term, resulting in no major distortion 
or scattering to smaller angles. 

The OWA of the image is set by the first Bessel term in the series in Eq. [2J We can 
get an approximate location of the OWA by noting that the first maximum of J n (x) is 
at approximately x = n. (A more precise bound may be found in IWatsonl (119441 ). which 
shows that the first maximum of J n occurs at n + 0.808618n 1/3 + C(ra~ 1/3 ).) The majority 
of the OWA peak comes from the first Bessel term in the perturbation series, which has 
magnitude: 

-R 



-2n 



< 



< 



< 



iXf 
4 

aj 

4 

a7 

2R 2 

a7 



I Jn {^f) ^^Irdr ) (2 cos (iV(0 - n/2))) 
( 2ixrp 



„ } "\—) rir 

r (2i\rp\ 



max 



o 

max 



J; 



N 



rdr 



Jn VW) 



(3) 



This bound is maximized when (27rrp) / (A/) ~ N, the first maximum. The maximum value 
of r is R = D/2, with D being the diameter of the pupil, so when p is at an angle of 
N/tt x X/D — or (Nf\)/(7rD) in distance units — then 



2nrp _ DN Xf 1 
~X~f ~ ni 2^lD\f 



N 



(4) 



5 



This provides a simple rule-of-thumb: 

N 

Powa ~ — (5) 

7T 

that is, the radius in A/D for the outer working angle with N holes is approximately N/ir; 
further iteration may be done locally 

One additional consideration is the interaction between the outer working angle 
structure and the aberrated light scattered outside the edge of the dark hole by the AO 
system. The electric field at the image plane is the convolution of the field from the 
aberrations and the PSF of the shaped pupil APLC; if the OWA is too small, then the 
convolution of a speckle in the halo around the dark hole and the APLC PSF will scatter 
light from the outer structure into the nominally dark region. A design goal is to have the 
OWA be greater than the longest distance across any part of the dark hole produced by the 
wavefront control system. If a system can correct up to ±/ m X/D in the x- and ^/-directions 
about the center of the PSF, where f m is the highest spatial frequency correctable by the 
AO system, then the longest distance within the dark hole is the diagonal between two 
opposing corners, with length 2y/2f m . Using the rule-of-thumb in Eq. [5j this means: 

N > 2nV2f m (6) 

that is, the pupil should have at least 2 r K\[2f m holes. For example, the GPI AO system can 
correct spatial frequencies up to 22X/D in a square about the core; the longest distance 
within the dark hole is ~ Q2.2X/D, which means the pupil should have at least 196 holes. 
(As with other design guidelines, this should be followed by a certain amount of local 
iteration to achieve the desired performance.) 

For a physical mask, simply introducing a number of holes is not good enough. 
Manufacturing requirements limit the size of any feature, so the design must be modified 
to take this into account. In addition, holes, and the bars between holes, should be larger 
than the wavelength of the light for which the mask is to be used, to ensure the equations 
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in this section remain valid and avoid breaking the assumptions we made in order to use 
scalar diffraction theory. 

To limit the minimum feature size, we "compress" the apodization by replacing A(r) 
with cA(r), where < c < 1. When this compressed apodization is converted into a binary 
mask, it creates a set of solid, wedge-shaped bars between the holes, ensuring a minimum 
width. Ideally, c should be as close to 1 as possible, as it decreases the throughput of the 
mask by decreasing the width of the holes in the mask. (See Fig. [3] for a visual depiction 
of the bars.) We also note that if these bars are sufficiently large, the pupil may be 
made as a free-standing pupil without a glass substrate. This would have the advantage 
of eliminating ghosting and chromatic effects from the mask, and potentially reducing 
aberrations introduced in the pupil. 

A second modification limits the size of the gaps between holes near the center, which 
can become very small for large numbers of holes. To do so, we truncate the hole sizes when 
they reach a certain minimum width. We then introduce a new set of holes in this region, 
with e.g. 50 holes instead of 200. We can determine how many holes are in the new region: 
let the truncation occur at a radius aR = aD/2, with < a < 1. (See Fig. [3] for a visual 
comparison of R and aR.) The perturbation expansion can be written in a similar form to 
Eq. m 

-^im, binary 

x(2cos (jN^-tt/2))) 

x(2cos(jiV 2 (0-7r/2))) (7) 
The second series will introduce perturbations around N 2 /tt x X/D, as shown above in Eq. 
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[3j The first series can be bounded similarly: 
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Here the maximum value of r is aR; we wish to choose a such that this Bessel function 
reaches its maximum at N 2 /iv x X/D, like the first series, so they both contribute at the 
same OWA location. This occurs at: 

2nrp 



iVi 



A/ 



-D A^2 Xf 1 

27ra 

2 n D Xf 



aN 2 



(9) 



As long as the number of holes Ni in the truncated region is greater than aN 2 , the truncation 
should not scatter additional light below the OWA. (Since these are approximations, some 
local iteration may be required.) Note that it is possible that the holes in the truncated 
region fall below the desired minimum size at some points; truncation can then be applied 
again, or multiple times, to create a new set of holes that are big enough. These truncation 
methods cannot be used for small troughs on the outer edge, however, since here the radius 
is that of the full aperture and the full number of holes is required to maintain the desired 
OWA. 



3. Results 

The design created for GPI is a 200-hole design, made with 50 holes in the truncated 
region. (See Fig. [2]for the continuous apodization and Fig. |4]for its 200/50-hole equivalent.) 
These parameters were chosen to maintain a 10/im minimum bar and hole size for a 12mm 
pupil over a 1.48-1.78/im band. We note that Eq. [2] and Eq. [3] are not dependent on the 
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size of the pupil; pupil size enters only in making modifications for manufacturing. GPI has 
a relatively small pupil, and so the hole truncation method is needed to keep the hole sizes 
up in the center of the star-shaped apodizer. 

To simulate the performance of the pupil designs under realistic conditions, we 
used a set of phase and amplitude aberrations that represent the static component of 
a typical ExAO system wavefront error. Dynamic aberrations were not included, since 
these result in a PSF halo whose speckle pattern averages out with time; we simulated 
the static aberrations that result in quasi-static speckles that limit ultimate sensitivity. 
The phase aberrations have an RMS of 5 nm from 0-3 cycles/pupil and 45 nm from 3-21 
cycles/pupil over the 7.8-m Gemini pupil, with a power law typical of manufactured optics; 
additional aberrations at subsequent optics are applied at 2.5-7 nm from 0-3 cycles/pupil 
and 0.5-21 nm from 3-21 cycles/pupil. Subsequently, near-perfect correction is assumed 
within the controllable range of spatial frequencies; residual phase offsets representing AO 
system calibration errors are reinjected at 5 nm from 0-3 cycles/pupil and 1 nm from 3-21 
cycles/pupil. 

The amplitude aberrations have an RMS of 0.1% — 0.3% percent at each optic. 
Amplitude aberrations are dominated by two effects: reflectivity variations of the primary 
mirror and Fresnel/Talbot effects, which cause a pure phase error on an intermediate- non- 
conj ugate optic to mix be tween phase and amplitude as the light propagates to the pupil. 



(See iMarois et al.l (j2008bl ) for further information.) 



The starshaped ap odizer is based on the current apodized design for GPI; see 



Macintosh et al 



(120081 ) for further details. The propagation between the incident plane and 



the Ly ot stop plane is done with the fast Lyot propagation algorithm created by Soummer 



et al. (ISoummer et al. 



20071 ); the PSF for the continuous apodized pupil is in Fig. [5l and 
the PSF for the starshaped apodizer is shown in Fig. [61 Radial contrast is calculated from 
the standard deviation of the PSF for all points in an annulus centered around a given 
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radial position. Spectral differencing may also be performed (see iMarois et al.l (120001 )). 
subtracting PSFs at multiple wavelengths and using the standard deviation calculation on 
those. In Fig. [7J we compare the performance of the nominal star-shaped design with the 
perfect apodized pupil and two additional starshaped apodizers without the truncation 
described above; the properties are summarized in Table [TJ Note that the 200/50-hole 
design performs comparably to the apodized and 200-hole designs, but with a minimum 
size on the holes and bars. The 50-hole design, while having relatively large holes, provides 
considerably poorer performance — this is a consequence of the interaction between the 
light beyond the AO frequency limit and the OWA, and emphasizes the need for a certain 
minimum number of holes. 



4. Conclusions and future directions 

We have demonstrated the theoretical performance of a shaped pupil in an APLC 
for the Gemini Planet Imager; the APLC performance with a starshaped apodizer is 
comparable to that with a perfect continuous apodizer, but with a minimum feature 
size of 10/im which is much larger than the wavelengths involved (typically 1.0 to 2.4 
microns). This technique thus avoids errors due to structures on the order of a wavelength. 
Random microdot pa tterns with dots larger than a wavelength are also a promising option 



Martinez et al 



2009|) 



The next step is to examine the performance of this design at the GPI Testbed at 
the American Museum of Natural History and also at the Laboratory for Adaptive Optics 
at University of California Santa Cruz. This pupil can be printed on glass or etched as a 
free-standing pupil; the latter may have the advantage of eliminating aberrations due to the 
glass substrate. We hope to confirm the theoretical performance both with and without the 
expected levels of aberrations due to propagation through turbulence. 
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Apodizer physical parameters 





Binary 


Binary 


Binary 


Perfect 




(200/50-hole) 


(200-hole) 


(50-hole) 


Apodization 


Diameter (mm) 


12 


12 


12 


12 


Obscuration ratio 


0.1432 


0.1432 


0.1432 


0.1432 


Mask size (A/D) 


5.6 


5.6 


5.6 


5.6 


Mask wavelength (/im) 


1.65 


1.65 


1.65 


1.65 


Number of main region holes 


200 


200 


50 


N/A 


Minimum feature size (holes) (/im) 


10.0 


3.5 


13.9 


N/A 


Minimum feature size (bars) (/im) 


10.2 


0.0 


0.0 


N/A 


Apodization compression factor 


0.90 


1 


1 


1 


Truncate troughs at inner edge? 


Yes 


No 


No 


N/A 


Number of holes in truncated region 


50 


N/A 


N/A 


N/A 



Table 1: Design parameters for the three binary designs and the continuous apodizer. 
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Fig. 1. — The parts of an APLC: 1. Entrance aperture with apodizer 2. Image plane 
hard-edged mask 3. Lyot stop 4- Final image plane 




Fig. 2. — The apodized pupil. 
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Fig. 3. — The binary starshaped pupil apodizer, with annotations. Top. The inner and 
outer regions, and their radii. The rectangle shows the area blown up below. Bottom. A 
close view of the binary apodizer, showing the bars and holes. 



11 




Fig. 4. — The binary starshaped pupil apodizer. The outer region has 200 holes, and the 
inner region has 50. 



Image plane slice after APLC at 1.5747um 

Image plane after APLC at 1 .5747um -4 i ■ ■ ■ n-> 1 1 — 
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Fig. 5. — The point spread function for the reference apodized pupil at 1.5747/xm. Aberra- 
tions are applied as described in Sec. Left. The entire PSF. Right. A slice through the 
center of the PSF. 
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Image plane slice after APLC at 1.5747um 

Image plane after APLC at 1 .5747um -4 i ■ ■ ■ n-> . . — 
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Fig. 6. — The point spread function for the starshaped binary pupil at 1.5747/xm. Aberra- 
tions are applied as described in Sec. [3] Left. The entire PSF. Right. A slice through the 
center of the PSF. 




Fig. 7. — The performance of three binary designs vs. the reference apodized pupil in the 
presence of typical aberrations. Left. Contrast after radial averaging. Right. Contrast after 
spectral differencing and radial averaging. 
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